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Abstract 

O 1 The Leutwyler-Smilga prediction regarding the (ir) relevance of the 

q ■ global topological charge for QCD in a finite box is subject to a test. 

£NJ . To this end the lattice version of a suitably chosen analogue (massive 

| 2-flavour Schwinger model) is analyzed in the small (ySm <C 1), inter- 

mediate (VT,m ~ 1) and large (FEm 3> 1) Leutwyler-Smilga regimes. 

■ The predictions for the small and large regimes are confirmed and il- 
[ — \ lustrated. New results about the role of the functional determinant in 

| all three regimes and about the sensitivity of physical observables on 

the topological charge in the intermediate regime are presented. 

<n : 

(N ■ 

o ■ 

o 1 Introduction 

o ■ 
o ■ 

■ One of the relevant concepts in an attempt to understand the mechanism 
of spontaneous breakdown of chiral symmetry^ in QCD is provided by in- 
stantons, i.e. topologically nontrivial solutions of the classical field equations 
which are both a local minimum of the classical action and localized in space- 
time. They are known to influence the local propagation properties of the 

; light flavours (Instanton Liquid Model 

^ ■ What is not so clear is whether the global concept of the number of 

instantons minus anti-instantons for QCD in a finite box 



plays a role, too. This is important, because in lattice studies involving 
dynamical fermions standard simulation algorithms tend to get "stuck" in a 
particular topological sector if the sea-quarks are taken sufficiently lightQ and 
obviously one would like to know whether this affects physical observables. 
In other words, the question is whether non-ergodicities of the sample w.r.t. 
the topological charge as visible e.g. from the r.h.s. of Fig. 1 tend to afflict 
measurements performed on such a sample. 

1 We shall only consider the theory with several light dynamical fermions {Nj > 2). 

2 The phenomenon was first observed for HMC with staggered quarks ||]. For Wilson 
type sea-quarks mobility between the topological sectors was reported to be sufficient for 
Mn/Mp > 0.56, but the same problem emerges if k is tuned sufficiently close to Kait H [§■ 
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Figure 1: LHS: Distribution of the unrenormalized field-theoretic topological 
charge v n& \ in an ensemble of 5000 SU (3) -configurations generated from <S\yiison 
at (3 = 6.1, after 6 cooling sweeps and with a 2-smeared charge operator (hgure 
taken from ^j). RHS: Time history, in units of MD time, of v na ,i for a HMC sim- 
ulation with a quadruple of staggered quarks at (3 = 5.35, m = 0.01; the analogous 
distribution lacks the symmetry v*-*—v and shows (u) ^0 (hgure taken from /fj/j. 

The correct sampling of the topological sectors is known to be crucial for 
quantities such as the rj' mass, which depend directly on the distribution of 
topological charges via the explicit breaking of the U(1)a symmetry ||. On 
the other hand, standard observables which do not relate to the U(1)a issue 
(M w , M K , M n , M p , V q q, . . .) are known not to depend on the global topological 
charge v if the four-volume V of the box is taken sufficiently large. An early 
investigation trying to assess the issue on a more quantitative level is the 
one by Leutwyler and Smilga |J . Their main observation is that for the case 
of pionic observables a complete answer can be given from purely analytic 
considerations for two extreme regimes of quark masses and box volumes, 
both of which involve the Leutwyler-Smilga (LS) parameter^ 

x = 1/Sm . (2) 

Their first statement is that in a sufficiently small box (x < 1) the question 
whether the algorithm managed to tunnel between the various topological 
sectors at a sufficient rate shouldn't emerge, because the partition function 
is completely dominated by the topologically trivial sector. In lattice lan- 
guage this amounts to the statement that the algorithm is not supposed to 
tunnel into a nontrivial sector anyways. Their second statement is that in 
the opposite regime (x ^> 1 plus a condition to be discussed below) the total 
topological charge is an irrelevant concept. This means, in lattice language, 
that the observable is unaffected by whichever "topological path" the simula- 
tion followed. From a lattice perspective, however, the analysis by Leutwyler 
and Smilga is not quite sufficient, as their paper doesn't state explicitly 
whether a simulation in the small x regime can be trusted, if the algorithm 
got stuck in a higher topological sector even though it was not supposed to 
get into it, and also because there is no statement about the situation in the 

3 V is the four- volume of the box, £ = — lim m ^o hniy_>oo (V'V') (this order) is the chiral 
condensate in the chiral limit and m is the (degenerate) sea-quark mass; note that both £ 
and m are scheme- and scale-dependent, but the combination is an RG- invariant quantity. 
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intermediate (x ~ 1) regime, which might be interesting to use in dynamical 
simulations. Last but not least the relevance of the second condition in the 
LS-analysis for large x (a condition which, as we shall see, is never fulfilled 
in QCD simulations) is not discussed. 

Below numerical results will be used to elucidate in which regimes of box 
volumes and quark masses a standard observable like the heavy quark po- 
tential is indeed insensitive to the global topological charge and by which 
mechanism a dependence sneaks in if the "insensitive" regime is left. I start 
with a brief review of the approach taken by Leutwyler and Smilga, followed 
by an argument why the massive multi-flavour Schwinger model is a suit- 
able testbed for numerical studies. Sections 4 through 6 present distinctive 
features by which the three LS-regimes differ from each other — both on 
a formal level (role being played by the functional determinant) and w.r.t. 
observables (Polyakov loop, heavy quark potential). The strategy used is 
to disentangle the complete sample into contributions from the individual 
topological sectors and to observe how the generating functional weights the 
fixed-z/ expectation values to form the physical observable [[/J. I conclude 
with a few remarks on the potential relevance of the results for future QCD 
simulations. 



2 Leutwyler- Smilga analysis for QCD 

Here I shall briefly summarize the considerations by Leutwyler and Smilga 
regarding the distribution of topological charges for QCD in a finite box ||. 

Leutwyler and Smilga start from the partition function for QCD with 
Nf>2 degenerate light (dynamical) flavours on an euclidean torus 

Z~J2 I det(p+m) N f e - (1 ^I GG (3) 

where the integration is over all gauge potentials in a given topological sec- 
tor and the quark fields have been integrated out. In (Q) the Vafa Witten 
representation [§] for the functional determinant 

det(p + mff = (m M li(X 2 +m 2 )) Nf (4) 

may be used which splits the latter into a factor representing the contribution 
from the zero modes of the Dirac operator p on the given background times 
a chain of factors ±iA+m from the paired non-zero modes. Upon using the 
nontrivial input for the first few non-zero eigenvalues 

a. = m (5) 

(which is just the familiar Banks Casher relation ||) the primed factor in 
the decomposition is seen not to depend on m if x<l, and hence 

Z v ~ m MN f {x < 1) . (6) 
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This means that in the chiral symmetry restoration regime (x < 1) topo- 
logically nontrivial sectors are heavily suppressed w.r.t. the trivial one — in 
particular for small quark masses (in lattice units) and large Nf. 

On the other hand, if the box volume is sufficiently large, the onset of 
chiral symmetry breaking gets manifest as the long-range Green's functions 
get dominated by the Goldstone excitations and the theory is effectively 
described by a chiral lagrangianfj 

Z~ / DUe~f LxpT , Lxpt = — (dUdlP) - -(UM^ + MU ] ) (7) 
Jsu(N f ) 4 2 

where the integration is over all smooth fields taking values in SU(Nf), and 
standard chiral perturbation theory (XPT) to order 0{p 2 ) has been usedQ. 
In the representation (Ut) U = exp(i\/2 4>/F) is the Goldstone manifold and 
M = diag(m, . . . , m) > the quark mass matrix; (. . .) means a flavour trace, 
and F is the pion decay constant to 0(p 2 ). The key observation by Leutwyler 
and Smilga is that if they take not only the box length sufficiently large for 
XPT to apply^ but also the quark masses sufficiently light 

1 « L «^r (8) 



Axpt ^ M v 

the collective (constant mode) excitations are no longer suppressed, but pro- 
vide the dominant contribution [II], i.e. ([?[) is approximately given by 



Z~ JdU e W(U Ml+MuS)/2 (Q) 

where the path integral has collapsed into a simple group integral given by the 
Haar measure. From (|9|) one gets the contribution of the topological sector v 
to the partition function in the standard way by substituting M -> M e w/N f 
and Fourier transforming from 9 to u, whence 

Z u ~ j\eJdU e-^e v ^ Uoe ~ iB/NfM ' +h - c ^ 2 

dU (detC/r e VE ^ Mt+h ' c -)/ 2 (10) 

U{N f ) 

where in the last line Uq and the factor e~ ie l N i have been combined into the 
U(Nf) matrix U . From (0) it follows that, if VEm is large, the dominant 
contribution to the group integral will come from the area where U is in 
the vicinity of the identity matrix, which means that the determinant is 
approximatively one and hence Z v is independent of v. A refined analysis 
shows that this is true for v 2 VEm/Nf, and the overall distribution is 

i? VY>m 
Z v ~ e ^> with (u 2 ) = -— (x > 1) . (11) 



4 References to all the subtle points (e.g. why in the particular case of the torus no 
boundary terms show up in Lxpt) are found in 

5 For a discussion in the framework of generalized chiral perturbation theory see Jl(| . 
6 Axpt denotes a QCD intrinsic low-energy scale, e.g. Aqcd, F^. 
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I shall attach a few comments: 

(i) It is worth noticing that the LS-classification of small, intermediate 
and large box volumes (i<l,i~l,a;>l) does not coincide with the usual 
classification on the lattice (L < M" 1 , L ~ M" 1 , L > M~ l ). Note also that 
one involves the masses of the sea-quarks, the other of the current-quarks. 

(ii) The analysis by Leutwyler and Smilga covers the symmetry restora- 
tion regime (k<1) and, on the other side, the regime where SSB is manifest, 
but the Goldstone bosons overlap the box (x^>l, L^CM^ 1 ). The latter con- 
dition — besides not being useful because it prevents^ extraction of pionic 
observables — is, in fact, completely inaccessible in present days simulations: 
Putting numbers into eqn. (§) one gets 

A XPT ~200MeV L~3fm M n ~ 20 MeV (12) 

which amounts to a sea-quark mass^ of the order of 70 keV. From this we 
conclude that, in order to be useful in dynamical QCD simulations, the result 
([TT| ) must turn out to rely exclusively on the condition i>1 being fulfilled, 
not on the r.h.s. of (^). In other words, our hope is that the latter turns out 
to be a purely technical condition, immaterial to the result (|TT|). 

(Hi) Comparing the announcement of the results by Leutwyler and Smilga 
in the introduction to their considerations as sketched above, one might 
worry because the latter aims at the relative weight of the topological sec- 
tors in the partition function, whereas the introduction mentioned possible 
^-dependencies of observables. From a formal point of view one may argue 
that this is not an issue, because the LS-analysis goes through if the theory 
is coupled to tiny external sources, and infinitesimal sources are sufficient to 
produce Green's functions and observables. Hence one expects that a strong 
dependence or approximate independence of the partition function w.r.t. v 
transfers into a sensitivity or immunity of correlation functions and observ- 
ables, but still one would desire to see this explicitly in numerical data. 

3 Adaptation to 2-flavour QED(2) 

The aim is to provide simulation data in the three LS-regimes i<1, x ~ 1 
and x^>l, while in addition fulfilling the usual condition L^>l/M n , i.e. the 
opposite of the r.h.s. of @. 

A toy theory where the LS-issue can be studied at greatly reduced costs 
(in terms of CPU power) while maintaining all essential features is the mas- 
sive multi-flavour Schwinger model, i.e. QED(2) with Nf>2 light degenerate 
fermions. The physics of this theory resembles in many aspects QCD (4) 
slightly above the chiral phase transition, i.e. the chiral condensate vanishes 
upon taking the chiral limit and the expectation value of the Polyakov loop is 

7 Unless the current-quark mass is taken considerably heavier than the sea-quark mass. 

8 Since M% ~ m, a sea-pion which is lighter than the physical pion by a factor 7 amounts 
to a sea-quark mass 49 times smaller than the phenomenological value of (m u + m ( j)/2. 
For quark masses the usual conventions (MS, /i = 2GeV) are adopted. 
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real and positive. Furthermore, the required formal analogy holds true: Like 
in QCD the gauge fields in the continuum version of the Schwinger model 



quantized on the torus fall into topologically distinct classes fll2| , i.e. they 
can be attributed a topological index. The latter agrees with the number of 
left- minus right-handed zero-modes of the massless Dirac operator, i.e. the 
index theorem holds true fl~2H . 

There is, however, an (apparent) problem: As we have seen in the previous 
section, the physics in the large LS-regime (i>1) is governed by the onset 
of SSB for the global SU(Nf)A flavour-group (though the box volume is still 
finite), but it is well known that in 2 dimensions no SSB with the associate 
production of Goldstone bosons takes place |T3| . 



The reason why this apparent deficiency proves immaterial relates to 
the fact that the multi-flavour Schwinger model (with Nf massless flavours) 
shows a second order phase transition with critical temperature T c = 0, a 
result due to Smilga and Verbaarschot . The theory is most appropriately 



seen as the limiting case of an extension where the axial flavour symmetry 
is explicitly broken, e.g. by tiny quark masses. In the strong-coupling limit[] 
m<^g the spectrum is found to involve a heavy particle with mass K 



M+ ~ v/iV>4= + °( m ) ( 13 ) 

V \/7T 



and N? — 1 light particles with mass |13], |16 



M_ ~ const ^/(AHD m N f /(N f+ i) (14) 



The state (13) is often called a "massive photon", but it is more appropriately 
seen as the analogue of the r/ in QCD, as it is the lightest flavour singlet and 
stays massive in the chiral limit. The states ([14]) are "Qwasi-Goldstones" 
and we shall call them "pions", but it is important to keep in mind that 
they differ from the usual "Pseudo-Goldstones" (as encountered in QCD) 
as they do not turn into true Goldstone bosons upon taking m — > 0; they 
rather become sterile if the chiral limitQ is performed. That conforms with 



Coleman's theorem |L3j which forbids the existence of massless interacting 
particles in 2 dimensions. The point I shall stress is this: As long as the 
axial flavour symmetry is explicitly broken — which, in the following, is 
true, as sea-quarks are taken massive — the "Quasi-Goldstones" are, for 
practical purposes, as good as "Pseudo-Goldstones" ; the difference shows 
only up upon taking the chiral limit. Hence even the large LS-regime, where 
long-range Green's functions are dominated by pionic excitations, finds its 
analogue in the (massive) multi-flavour Schwinger model. 

One more practical issue has to be resolved: The original definition (0) 
of the LS-parameter x involves the quantity S, i.e. (the absolute value of) 

9 Note that in 2 dimensions the charge has the dimension of a mass. 
10 Analogous phenomena are observed, if the axial flavour symmetry is broken by bound- 
ary conditions rather than quark masses: The chiral condensate (i.e. the "would be order 
parameter") vanishes with a power-law dependence as the spatial box length is sent to 
infinity at zero temperature, but exponentially fast if the temperature is non-zero |l7|. 
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the chiral condensate in the chiral limit. The latter, however, is exactly zero, 
due to Coleman's theorem. A practical way out of this is based on Smilga's 



observation that upon combining (|14|) with the expression [16 



| I ^ const £ 2 /W +1 ) m W-D/(AHD (15) 



for the chiral condensate one gets the relation |L6| 

m\(M)\ ~ , 0,388 - Ml (16) 
(2.008) 2 77 y 1 

where I have already plugged in the nonuniversal (i.e. ^-dependent) con- 
stants for the 2-flavour case fL8| . Relation ( |16|) being the QED(2) analogue 



of the Gell-Mann-Oakes-Renner PCAC-relation 

mS ~ \MlFl (VN f ) (17) 

in QCD means that I may define the LS-parameter in the lattice studies 
presented below as 

* s msr v * (N '- 2) (18) 

as one may rewrite it in the case of QCD (to leading order in XPT) 

x= X -VMlFl (ViV/) (19) 

in a way which involves just physical degrees of freedom^ 

In order to investigate the LS-issue, I have chosen to implement the 
Schwinger model with a pair of (dynamical) staggered fermions, using the 
Wilson gauge action S^auge = P — cos6> n ). Since the staggered Dirac op- 
erator is represented by relatively small matrices (see below), I have decided 
to compute the determinants exactly, using the routines ZGEFA and ZGEDI 
from the LINPACK package. 

I compare the three regimes x-C l,x ~ l,x ^> 1 to each other using 
three dedicated simulations: I have chosen to vary the box volume at fixed 
(3 = l/(ag) 2 = 3.4 and fixed staggered quark mass m = 0.09 (everything 
in lattice units). The three regimes are represented by the three volumes 
V = 8 x 4, V = 18 x 6, V = 40 x 10, where always periodic boundary conditions 
in the first (spatial) direction and thermal boundary conditions in the second 
(euclidean timelike) direction have been used. The LS-parameter ( jT8|) takes 
the values x ~ 0.33, x ~ 1.12, x ~ 4.16 respectively, while the pion (pseudo- 
scalar iso-triplet) has a (common) massQM 7r = 2.008 • 0.542 1 / 3 0.09 2/3 = 0.329 
and hence a correlation length £„- = 3.04 as to fit into the box (almost) in all 

n Note that in order to determine a; in a lattice study through ( |l9| ) one has to plug 
in Mir and F„ of the sea-pion, i.e. in practice of a pion constructed from current-quarks 
having exactly the same mass as their sea-partners. 

12 The formula used is the prediction (|lj), which seems adequate since m/g~0.166-C 1. 
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Figure 2: Distribution of (left) and f geo (right) in the small, intermediate 
and "large" Leutwyler-Smilga regimes, respectively. 



three regimes - as is usual in a lattice context. Each time 3200 decorrelated 
(w.r.t. the plaquette) configurations have been generated. 

For the type of investigation I am aiming at configurations must be 
assigned an index v. I have implemented both the geometric definition 



j- log Un and the field theoretic definition z/f t h = k z/ n 



: E sin 9 n 



with the renormalization factor k~1/(1 — {S gange ) / j3V) |19| . A configuration 
is assigned an index only if the geometric and the field-theoretic definition, 
after rounding to the nearest integer, agree. This turned out to be the case 
for 3197, 3163, 2818 configurations, i.e. on a 99.9%, 98.8%, 88.1% basis on 
the small/intermediate/large lattice. This fraction being so high means that 
in practice an assignment can be done without cooling for the overwhelm- 
ing majority of configurations. The remaining ones are just not assigned an 
index; they are used to compute the unseparated observable (e.g. the "com- 
plete" heavy quark potential; see below), the latter, however, turns out to 
be unaffected by whether they are included or not. 

Having simulation data which are supposed to reflect the situation in the 
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small/intermediate/large LS-regimes respectively, one may want to check the 
overall distribution of topological charges as encountered in the three runs. 
As one can see from the r.h.s. of Fig. 2, the overall distribution is noticeably 
different in the three regimes. The predictions by Leutwyler and Smilga for 
the two extreme cases (x <C 1 and x ^> 1) seem to be well reproduced: For 
small x the partition function Z (and hence the sample) is dominated by 
Zq, the contribution from the topologically trivial sector, whereas for large 
x the distribution gets broad and seems compatible with a gaussian. It is 
worth mentioning that the latter holds true even though the pion did not 
overlap the box; rather the usual relation applied. This is a first 

indication that the r.h.s. of (|8|) in the LS-analysis might indeed be a purely 
technical condition, immaterial to what they claim, namely that in the large 
x regime the topological charge is "an irrelevant concept" or, in other words, 
that physical observables do not depend on v if x ^> 1 || . 

4 Sectoral reweightings and renormalizations 

Since the LS-issue is peculiar to the full (unquenched) theory, an attempt 
to understand by which mechanism the three regimes differ from each other 
may lead one to investigate how the functional determinant or specifically 
its contribution to the total action per continuum-flavour 

Sfermion = - log(det (p + m) ) + const (20) 

relates to the contribution from the gauge field and, in addition, how such a 
relationship might depend on the topological charge of the background. The 
idea is thus to study such a relationship sectorally, i.e. after the complete 
sample has been separated into subsamples with a fixed topological charge. 
Since the theory (at 6 = 0) is P/CP-symmetric, the configurations from the 
sectors ±z/ may be combined into a single subsample characterized by \v\. 

As one can see from the rightmost column of Fig. 3, the contribution per 
continuum flavour to the total action, Sf er mion, shows a positive correlation 
with the gauge action .Sgauge- This means that, on a qualitative level, switch- 
ing on the functional determinant amounts to an increase of (3, which is a 
well-known feature in (full) QCD too. However, the correlation is weak, i.e. 
the "cigars" in the unseparated plots are thick. A key observation is that the 
correlation improves, if one separates the sample into subsamples with fixed 
as is done in Fig. 3. It is worth noticing that in general (i.e. without 
specifying the LS-regime) the best linear fit for Sfermion as a function of S gmge 
depends (in offset and slope) on \v\. This means that, to a first approxima- 
tion, the functional determinant results in an overall suppression of higher 
topological sectors w.r.t lower ones and to a sectorally different renormaliza- 
tion of (3. As we shall see, these two features constitute the place where the 
three LS-regimes differ from each other. 

First, we shall consider the offset of the regression line in the v = ±1 
sectors compared to the line in the sector with z/ = 0, as this is a measure how 
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much extra suppression the first topologically nontrivial sector receives on 
average compared to the trivial one due to the functional determinant. From 
the first and second column in Fig. 3 one sees that this offset (and in general 
the offset between neighboring sectors) decreases with increasing x. This 
provides an explanation for the broadening of the ^-distribution inherent 
to passing from the small to the intermediate or from the intermediate to the 
large LS-regime (cf. Fig. 2). For x < 1 the offset (5f ermioil )i - (S icriaion ) is 
larger than the fluctuations in Sfermion within the zero-charge sector, and this 
means means that in the small LS-regime the functional determinant acts 
as an almost strict constraint to the topologically trivial sector. For x ~ 1 
the offset has roughly the same size as typical fluctuations of Sfermion in the 
sector with v = 0; as a result of this one can say that in the intermediate 
regime the functional determinant acts as a "soft constraint" to low-lying \u\ 
values. For x ^> 1 finally, (Sfermion) i — (Sfc rm ion)o is small compared to the 
fluctuations in each one of these sectors; hence no such effective description 
proves appropriate in the large LS-regime. 

Passing on to the slope of the regression line in the correlation (Sfermion) 
versus (S gauge ), we notice from Fig. 3 that this slope decreases (in general) 
with increasing \v\. This means that the effective renormalization 

/3<?) ptg = (1 + Nf sloped)) (21) 

due to the functional determinant is (in general) stronger for sectors with low 
absolute value of the net topological charge. This effect of sectoral depen- 
dence of the effective renormalization factor for (3 diminishes with increasing 
x; eventually, in the large x regime, it is completely gone. 

Thus, the large LS-regime is unique, as the determinant results only in 
an overall suppression of higher topological sectors, but not in a sectoral 
dependence of f3 e s (and hence, in 4 dimensions, of the lattice spacing in 
physical units). This seems consistent with the claim by Leutwyler and 
Smilga that physical observables do not depend on v in the large x regime, 
but it is of course not a sufficient condition for the latter to be true. 

5 Sectoral heavy quark free energies 

A quantity which is easily accessible in lattice studies and which might help 
to elucidate the physical meaning of the LS-classification is the Polyakov 
loop, i.e. the trace of a chain of link variables which winds once around the 
torus in the euclidean timelike direction. Its logarithm represents (up to a 
factor) the free energy of an external ( "heavy" ) quark which is brought into 
the system without possibility to influence it through back-reaction. As in 
the previous section the analysis shall be done sectorally, i.e. on classes of 
configurations with a fixed value of \v\. 

The result of such an analysis for the Polyakov loop is displayed in Fig. 4. 
As one can see from the rightmost column, the expectation value (L) on the 
"complete" (unseparated) sample decreases if the box volume is increased. 
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What is more interesting is the following: The physical (unseparated) expec- 
tation value may be understood as a weighted mean 

(L)= J2 PH(L) H (22) 

\u\>0 

of sectoral expectation values (L)\ v \, where the factor p\ v \ reflects the com- 
bined weight of the sectors ±v in the histogram on the r.h.s. of Fig. 2. From 
the first line in Fig. 3 one notices that the sectoral expectation values for 
\v\ =0,1 in the small LS-regime differ quite drastically on the scale set 
by the physical expectation value, i.e. |(L)|o|,x<i — (-^}|i|,x<i|/(^)cc<i — 0.7. 
In the intermediate regime the analogous ratios for neighboring sectors (i.e. 

\{ty\o\,x~i-{ty\i\,x~i\/{L)x~i and | {L)\ 1 \ )X ~i-{L)\ 2 \,x~i\/ (L) x ~x) happen to be 
smaller, but comparing the highest to the lowest available \u\, one finds again 
\(L)\o\,x~i- (L>)\2\,a;~i\/(L)x~i - 0.7. Finally, in the large LS-regime, there 
is an almost perfect consistency between sectoral Polyakov loop expectation 
values for low-lying \u\ (i.e. \u\ G {0, 1, 2}), but there is no numerical evidence 
in favour of an overall consistency; if one builds the analogous ratio for the 
lowest and highest \u\ available, the quantity |(L)|o|,a;>i — (■^)|3|,x»i|/(-^)x»i 
is found to be again of order one. 

The Polyakov loop is the first example of a physical observable for which 
we see a difference between expectation values in neighboring sectors (for 
sufficiently low \u\ values) disappear upon taking the large x limit. However, 
even in the large LS-regime no supporting evidence has been found that (L) v 
might be completely independent of v. It seems that (L)\ v \ proves sensitive 
on topology even in the large LS-regime when jumping from |^| m j n to |^| max - 



6 Sectoral heavy quark potentials 

The second observable which has been evaluated for each sector separately is 
the heavy quark potential as determined from the correlator of two Polyakov 
loops — see Fig. 5 and Tabs. 1-3. The first thing to notice is that in each 
regime the physical potential (the one in the rightmost column in Fig. 5) 
bends to the right. This curvature is a direct indication of the screening 
brought by the dynamical fermions; the quenched potential is strictly linear^. 

In the small LS-regime the sectoral potentials for u — and u = ±l differ 
quite drastically. In this case the potential in the topologically trivial sector 
agrees (on a 1 a basis, cf. Tab. 1) with the potential on the full sample — as 
expected, since for x <C 1 the partition function is dominated by Z . What 
might come as a surprise is that the error-bars in the first plot in the top 
line of Fig. 5 are smaller than those in the rightmost one, even though the 
configurations analyzed in the former case represent a subset of those used in 
the latter case. The reason, as we have seen, is that the complementary subset 
(the one comprising the configurations with z/ = ±l) provides measurement 

13 Unlikc in QCD(4) there is no 1/r short-distance tail; the one-photon exchange con- 
tributes a linear piece in 2 dimensions. 
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Table 1: Sectoral heavy quark potentials and unseparated (physical) poten- 
tial at spatial separations 1..3 in the small Leutwyler Smilga regime (x<^l). 
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Table 2: Sectoral heavy quark potentials and unseparated (physical) poten- 
tial at separations 1..4 in the intermediate Leutwyler Smilga regime (x~l). 
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Table 3: Sectoral heavy quark potentials and unseparated (physical) poten- 
tial at spatial separations 1..5 in the large Leutwyler Smilga regime (x^>l). 



data which are inconsistent with the data won from the topologically trivial 
subset, but has a total weight which is too small to affect the complete 
ensemble average significantly (cf. the first plot in Fig. 2). Hence, in the 
small LS-regime, the tiny contribution from the higher topological sectors 
primarily adds some noise to a standard observable like the heavy quark 
potential. The practical recipe to get rid of this effect is to do by hand (for 
sufficiently small x) what the functional determinant does anyways in the 
limit x— *0: Cut away the higher topological sectors ! 

In the intermediate LS-regime the difference between the two sectoral 
potentials for z/ = and z/ = ±l is smaller than in the previous regime, but 
the sectoral potentials are typically several a away from each other and even 
the potential in the topologically trivial sector lies significantly below the 
unseparated physical potential (cf. Tab. 2). From this we see that the sector 
with z/ = does not always yield the right result. The intermediate LS-regime 
shares thus with the small x regime the uneasy property that different sectors 
give inconsistent values for physical observables, but unlike in the small x 
regime, several sectors yield sizable contributions to the complete (physical) 
expectation value — hence perfect ergodicity w.r.t. v is mandatory. 
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In the large LS-regime differences between expectation values measured 
on subsets of definite topological charge diminish further, if the latter differs 
only by one ("neighboring sectors"). On the other hand, in the large x regime 
there is a multitude of topological sectors yielding sizable contributions to 
the physical expectation value, hence an obvious question is which effect will 
win in the limit x^oo. Hence, what we would like to know is whether the 
difference between the heavy quark potential on the subsets with the lowest 
and the highest (reasonably populated) \v\ will fade away, stay constant or 
grow if x increases unlimitedly. From Tab. 3, one gets the impression that the 
difference between the two extreme fixed- \v\ expectation values persists in the 
large LS-regime. It seems that the statement by Leutwyler and Smilga that in 
the large x regime the topological charge has vanishing influence on physical 
measurements is reproduced for sufficiently low \u\ values (in our case e.g. 
\v\ e {0, 1, 2} but not necessarily for \v\ =3). It is also worth noticing that the 
topologically trivial potential lies at least 0.5 a below the physical potential. 
Hence, a constraint to the topologically trivial sector seems permissible for 
x >> 1, though in our simulation a constraint to the sectors \v\ = 1 would 
have been a better choice. Of course it is difficult to get definitive answers 
from numerical data; it may well be that our LS-value (x ~ 4.16) is still 
too small for the large x behaviour to be fully pronounced. Nevertheless, 
the situation seems to be analogous to what we have seen in the case of 
the heavy quark free energy: Any differences between expectation values of 
neighboring sectors (with sufficiently low \u\) disappear upon taking the large 
x limit. However, even in the large LS-regime no supporting evidence has 
been found that {V q q) v might be completely independent of v. 

7 Summary 

The massive multi-flavour Schwinger model has been used to test the state- 
ments by Leutwyler and Smilga about the (ir-)relevance of the topological 
charge for QCD in a finite box. It is expected that an analogous study for full 
QCD will result in figures which look, as far as qualitative issues go, exactly 
the same as the ones presented here and that, for this reason, the statements 
below are valid for QCD too. The key points are the following: 

(i) The three LS-regimes show distinctive properties already on a formal 
level, by the way how the contribution per flavour to the total action (the 
sign-flipped logarithm of the one-flavour functional determinant) correlates 
with the contribution from the gauge field. Introducing the fermion deter- 
minant at fixed (3 results in every regime in an overall suppression of higher 
topological sectors w.r.t. low-lying ones. For i C 1 this suppression is so 
strong that the functional determinant effectively acts as a constraint to the 
topologically trivial sector. In addition, the sea-flavours effectively result in 
a sectoral multiplicative renormalization of (5 with a factor greater than 1. 
For small and intermediate x this factor is found to decrease as a function 
of while in the large LS-regime the renormalization factor seems to be 
uniform for all topological sectors. 
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(ii) For the overall distribution of topological charges the predictions by 
Leutwyler and Smilga are well observed: In the regime x <C 1 the partition 
function (and hence the sample) is entirely dominated by the topologically 
trivial sector, i.e. the histogram of topological charges essentially consists of 
a delta-peak at v — 0. In the intermediate regime (x ~ 1) the distribution 
is narrow but nontrivial, i.e. sizable contributions stem essentially from the 
sectors u — 0, ±1, ±2. In the regime i>1 the topological charge distribution 
gets broad and is approximated by a gaussian with width (z/ 2 ) ~ x/Nf. 

(Hi) In the small LS-regime (x <C 1) the sectoral expectation value of an 
observable in the first topologically nontrivial sector (typically won from very 
few configurations) differs drastically from the analogous expectation value 
in the topologically trivial sector (represented by the overwhelming majority 
of configurations). This means that a simulation getting stuck at z/ = ±l (or 
even higher) is absolutely disastrous to the result, i.e. tiny non-ergodicities 
w.r.t. the topological charge render an overall sample average meaningless. 
The good news is that in this regime the correct distribution is known and 
hence there is a simple remedy: A strict constraint to the topologically trivial 
sector proves beneficial. 

(iv) The intermediate LS-regime (x~l) is characterized by the distinctive 
feature that several mutually inconsistent sectors yield sizable contributions 
to a standard (i.e. not relating to the U(1)a issue) physical observable. This 
means that the intermediate regime is the one which is most delicate to 
simulate in: There is a strong sectoral dependence of physical observables, yet 
there is no simple recipe how to guarantee the correct sectoral weighting. In 
other words: Physical measurements on a sample won in the intermediate LS- 
regime do rely on the simulation algorithm having achieved perfect ergodicity 
w.r.t. the topological charge. 

(v) In the large LS-regime (i>1) sectoral averages seem consistent be- 
tween neighboring topological sectors (for sufficiently low \v\), but no evidence 
has been found that measurements might be completely independent of v. It 
seems that jumping (within a well-distributed sample) from the lowest to the 
highest accessible \v\ may affect sectoral averages even in the large x regime. 
This, if correct, means that the statement by Leutwyler and Smilga accord- 
ing to which "the topological charge is an irrelevant concept in the large x 
regime" || should to be interpreted in the following way: In the large LS- 
regime physical observables are immune to modifications of the v-histogram 
which are small compared to its natural width (x/Nf) 1 / 2 ; they may, however, 
be sensitive to distortions which go beyond this limit. From a lattice per- 
spective the implication is that the result of a full QCD simulation in the 
large x regime may be trusted without hesitation as long as the algorithm 
did not get stuck in a sector with unreasonably high say at u = 10 if the 
simulation was in the regime x~40 and Nf = 2. 

(vi) In this work, the transition from the small to the intermediate and 
the large LS-regime has been made by adding more sites to the grid while 
keeping (3, m fixed. It would be interesting to test the universality property, 
i.e. the claim that it is only the product VTjm || which decides on the regime. 
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(vii) It is worth emphasizing that in the large x simulation used in this 
work the predictions by Leutwyler and Smilga turned out to be fulfilled even 
though the condition in the r.h.s. of (§) according to which the pion would 
overlap the box was not obeyed, but rather the opposite relation l/M n <^L 
would hold true, i.e. the pion would fit into the box, as is usual in a lattice 
context. Hence the r.h.s. of (H) seems to be a purely technical condition, 
immaterial to the result of the LS-analysis in the large x regime. 

(viii) From the fact that the massive multi-flavour Schwinger model fol- 
lows the LS-classification even though it resembles QCD slightly above the 
chiral phase transition one is led to suspect that this classification might 
be more general, even within QCD, than its original derivation indicates. In 
other words: The conjecture is that the LS-classification may prove useful not 
only in the broken phase but also in the high-temperature phase of QCD - 
as long as one stays in the immediate vicinity of the phase transition, i.e. as 
long as pions keep being visible (i.e. keep dominating the long-range Green's 



functions with pion-type quantum numbers; for a recent study see e.g. [21]). 
This, if correct, would imply that the version of the LS-parameter as given 
in eqn. flI5|) is more general^ then the original version (0) . 

(ix) The results of this investigation may be seen as an additional a- 
posteriori justification of the resources phenomenological lattice groups have 
been provided with in order to be able to study (full) QCD in the large x 
regime[5 In addition, they support attempts to develop algorithms which 



decorrelate the topological charge faster than standard HMC (see e.g. j22 



as this is necessary for trustworthy simulations in the intermediate LS-regime, 
towards which one moves when comparing different runs with increasing K sea 
in a fixed physical volume, or even more so, if (3 is kept fixed. The fact that 
even in the large x regime full QCD simulations may yield unreliable results 
upon generating extreme distortions in the topological charge distribution 
should be taken as a motivation to monitor v in all phenomenological studies 
(i.e. not just with an observable relating to the XJ{V)a issue) and, for the time 
being, as a warning against pushing the sea-quark mass too light^ (even when 
increasing the volume together with K sea so as to keep the product x = VEm 
fixed and large compared to 1). The latter constraint is certainly not a 



disaster, since there is a well-elaborate framework [E3l] which allows one to 



extrapolate from samples generated with a somewhat larger sea-quark mass 
down to the physical M n /M p . Hence it seems that our results do represent 
a little caveat but no real obstacle to future progress in Lattice QCD. 
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